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ABSTRACT 

We present and analyze Far Ultraviolet Spectroscopic Explorer (FUSE) observations of six solar analogs. These 
are single, main-sequence GO-5 stars selected as proxies for the Sun at several stages of its main-sequence lifetime 
from ~130 Myr up to ~9 Gyr. The emission features in the FUSE 920-1180 A wavelength range allow for a 
critical probe of the hot plasma over three decades in temperature: ~ 10 4 K for the H I Lyman series to ~ 6 • 10 6 
K for the coronal Fe XVIII A975 line. Using the flux ratio C III A1176/A977 as diagnostics, we investigate the 
dependence of the electron pressure of the transition region as a function of the rotation period, age and magnetic 
activity. The results from these solar proxies indicate that the electron pressure of the stellar ^10 5 -K plasma 
decreases by a factor of ^70 between the young, fast-rotating (P m = 2.7 d) magnetically active star and the 
old, slow-rotating (P rot ~ 35 d) inactive star. Also, we study the variations in the total surface flux for specific 
emission features that trace the hot gas in the stellar chromosphere (C II), transition region (C III, O VI), and 
corona (Fe XVIII). The observations indicate that the average surface fluxes of the analyzed emission features 
strongly decrease with increasing stellar age and longer rotation period. The emission flux evolution with age or 
rotation period is well fitted by power laws, which become steeper from cooler chromospheric (~ 10 4 K) to hotter 
coronal (~ 10 7 K) plasma. The relationship for the integrated (920-1180 A) FUSE flux indicates that the solar 
far-ultraviolet (FUV) emissions were about twice the present value 2.5 Gyr ago and about 4 times the present 
value 3.5 Gyr ago. Note also that the FUSE/FUV flux of the Zero-Age Main Sequence Sun could have been 
higher by as much as 50 times. Our analysis suggests that the strong FUV emissions of the young Sun may have 
played a crucial role in the developing planetary system, in particular through the photoionization, photochemical 
evolution and possible erosion of the planetary atmospheres. Some examples of the effects of the early Sun's 
enhanced FUV irradiance on the atmospheres of Earth and Mars are also discussed. 

Subject headings: stars: late-type — stars: atmospheres — stars: chromospheres — stars: coronae — stars: 
activity — ultraviolet: stars 



1. INTRODUCTION 

The Sun's magnetic activity is expected to have greatly de- 
creased with time (Skumanich 1972; Simon, Boesgaard, & Her- 
big 1985; Dorren & Guinan 1994; Guinan, Ribas, & Harper 
2002) as the solar rotation slows down because of angular mo- 
mentum loss in the stellar wind and resultant reduction in mag- 
netic dynamo-related activity. The study of the young Sun's far- 
ultraviolet (FUV) fluxes using solar proxies provides important 
diagnostics for the state of the younger Solar System and the 
physics of the much more active early Sun. The comprehensive 
"Sun in Time" project, begun in 1988 (Dorren & Guinan 1994), 
focuses on the study of the long-term evolution of the outer at- 
mosphere of an early G star, from the zero-age main sequence 
to the terminal-age main sequence. A crucial component of the 
program is a carefully-selected sample of nearby solar analogs 
(Giidel, Guinan, & Skinner 1997) with spectral types confined 
between G0-5 V and with well-determined physical properties 
(including temperatures, luminosities, and metal abundances). 
We have obtained extensive photometry (UBVRI) of these stars 
from ground-based observatories for over a decade to determine 
their rotation periods, investigate starspots and possible activity 
cycles. In addition, we have been able to estimate the stellar 



ages by making use of their memberships in clusters and mov- 
ing groups, rotation period-age relationships, and, for the older 
stars, fits to stellar evolution models. 

The sample of solar proxies within the "Sun in Time" pro- 
gram (see Guinan, Ribas & Harper 2002) contains stars that 
have masses close to 1 Mq and cover most of the Sun's main 
sequence lifetime from ~130 Myr to ^9 Gyr. Basically, the 
program stars have similar convective-zone depths and chemi- 
cal abundances to those of the Sun and vary only by their age 
and rotation periods (P, ot ), and hence dynamo-generated mag- 
netic activity, i.e., chromospheric, transition region (TR), and 
coronal emissions. In essence, we use these one solar mass 
stars with P 10t between 2.7 and ~35 days as laboratories to study 
and test solar and stellar dynamo theories by varying essentially 
only one parameter: rotation period. 

The "Sun in Time" is a comprehensive and multi-frequency 
program that addresses a variety of topics: study of short- 
and long-term magnetic evolution; physics and energy trans- 
fer mechanisms of the chromosphere, transition region (TR) 
and corona; evolution of the spectral irradiances of the Sun 
and their effects on paleo-planetary environments and atmo- 
spheres. To this end we utilize observational data spanning al- 
most the entire electromagnetic spectrum, which allows us to 
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probe the structure of stellar/solar atmospheres. The FUV has 
not been readily accessible until the successful launch of the 
Far Ultraviolet Spectroscopic Explorer (FUSE). FUSE allows 
a critical probe of the hot plasmas over nearly three decades 
in temperature: e.g., ~ 10 4 K for the H I Lyman series (Ly/3, 
Ly 7 , ...) through O VI AA1032, 1038 at - 3 • 10 5 K, and the 
recently identified Fe XVIII A975 coronal emission feature at 
~ 6 • 10 6 K. Several of the strongest emission features, such as 
C III AA977,1 176 and O VI AA1032,1038, originate in TR plas- 
mas, and are pivotal for understanding the mechanisms of chro- 
mospheric and coronal heating. Thus, spectrophotometry with 
FUSE fills the energy gap between IUE/HST (1160-3200 A) 
and EUVE (80-720 A) to yield a complete and comprehensive 
picture of a solar-type star's atmosphere at different ages and 
rotation periods. 

In this paper we present and analyze FUSE observations of 
six bright solar analogs that span a wide range of rotation pe- 
riods (and ages) from 2.7 days (^130 Myr) to ^35 days (<~9 
Gyr). The stars, whose properties are listed in Table 1 along 
with the Sun's, are: EK Dra, ir 1 UMa, k 1 Cet, (3 Com, (3 Hyi, 
and 16 Cyg A. The present study focuses on two particular as- 
pects of the "Sun in Time" program that can be effectively ad- 
dressed with FUV observations: 1.) The flux ratios between 
the C III A 1176 and C III A977 lines (both within the FUSE 
spectral range) yield empirical measures and estimates of the 
electron pressure (P e ) of the TR; 2.) Integrated fluxes for se- 
lected emission lines can be used to study the evolution of the 
solar FUV irradiances at different ages and rotation periods. 

2. OBSERVATIONS 

FUV observations of the six targets discussed here were ob- 
tained with FUSE. The FUSE instrument consists of four spec- 
troscopic channels (LiFl, LiF2, SiCl, and SiC2) that cover a 
combined wavelength range 905-1 187 A. Detailed information 
on the FUSE mission and its on-orbit performance are provided 
by Moos et al. (2000) and Sahnow et al. (2000). 

The observations of the targets were secured during FUSE 
Cycles 1, 2 and 3. Spectra of three stars (k 1 Cet, f3 Com, and 
(3 Hyi) were taken as part of the Guest Investigator program 
A083 (RI. Guinan), one star (tt 1 UMa) was observed as part 
of the program B078 (RI. Guinan), and to additional stars (EK 
Dra and 16 Cyg A) were acquired within program C102 (RI. 
Guinan). The large aperture (LWRS) was employed in all cases 
to minimize photon losses due to telescope misalignments, and 
this caused a fairly severe contamination of day-time spectra 
from geocoronal emission lines and scattered solar radiation. 
Information on the dates and exposure times of the observations 
is provided in Table 2. 

The FUSE data were retrieved from the Multimission 
Archive at Space Telescope (MAST), and then re-calibrated 
with CALFUSE 2.0.5 4 . For each star, the time-tagged data for 
the sub-exposures were combined into a single dataset. The 
spectral extraction windows were individually aligned to ensure 
correct extraction of the stellar spectra. We employed the CAL- 
FUSE 2.0.5 event burst screening and default background sub- 
traction. The data were then calibrated and screened separately 
for orbital night and day, allowing us to estimate any contami- 
nation of the stellar spectra from "airglow" and scattered solar 
light. The scattered light was found to be negligible compared 
to the stellar source, except for the faint FUV target 16 Cyg A, 
which could suffer from some contamination especially in the 

4 http://fuse.pha.jhu.edu/analysis/calfuse.html 



O VI features. 

The calibrated spectra are oversampled so they were re- 
sampled by binning over four wavelength intervals, thereby in- 
creasing the signal-to-noise ratio. A comparison of the emis- 
sion line fluxes from different channels, e.g., C III 977 A SiC2A 
and SiClB, provides an independent check on data artifacts, 
and whether the star was in both apertures for the same length 
of time. 

3. ANALYSIS 

For illustration we show the extracted spectrum of k 1 Cet 
in Figure 1. Note that we only show the night-time data that 
presents less contamination by terrestrial airglow and geocoro- 
nal emission. Thus, all features in Figure 1 are of stellar ori- 
gin, although the H I Lyman series still has some geocoronal 
contamination. The FUV spectra of the other four targets are 
qualitatively similar, and only the spectrum of EK Dra and tt 1 
UMa are somewhat noisier. For these two stars the increased 
noise arises from the failure of the SiC2A and LiF2A observa- 
tions and the consequent use of the less sensitive segment B of 
detector 1 (SiClB and LiFIB) for those parts of the spectra. As 
can be seen in Figure 1, the continuum level is negligible and 
the flux is in the form of emission lines. Individual channels 
are represented in the figure because the FUSE effective area is 
heavily dominated by only one channel at each wavelength. In 
this case, SiC2A, LiFIA, and LiF2A spectra are shown in the 
upper, central, and lower panels, respectively. 

Studies of FUSE observations of late-type stars, AB Dor and 
Capella, have been presented by Ake et al. (2000) and Young 
et al. (2001), respectively. The authors carried out thorough in- 
vestigations of the chromospheric and TR emission features in 
the FUSE wavelength range and provided identification spec- 
tra. A comparison with our sample spectrum in Figure 1 indi- 
cates a striking resemblance to these spectra, except for small 
to moderate changes in the relative strengths of some features. 
We have thus been able to identify the most significant fea- 
tures in our spectra by using the results of Ake et al. (2000) 
and Young et al. (2001). Among the strongest FUV emission 
lines in coronal spectra of late-type stars are the C III singlet at 
977.020 A, the C III multiplet around 1 176 A, and the O VI dou- 
blet at 103 1 .925 A and 1037.614 A. We focus here primarily on 
the analysis of stellar features that are not significantly affected 
by geocoronal contamination and day-time spectra can be used 
without any loss of accuracy. We thus combined the night-time 
and day-time data by using a exposure-time weighting scheme. 

All of our FUSE spectra show clean profiles for these lines 
that permit reliable integrated flux measurements, which were 
carried out by means of NOAO/IRAF tasks. Different meth- 
ods (plain addition, Gaussian fitting) were employed to measure 
the integrated line fluxes, study the associated uncertainties, to 
check the line profiles, and also to assess the possible presence 
of interstellar absorption (see below for further discussion). The 
total integrated fluxes for each star and their intrinsic uncertain- 
ties (random errors) are provided in Table 2. 

The ratio of the fluxes for the C III A1176 and A977 emis- 
sion features is provided in Table 6. This ratio, which yields 
an excellent electron pressure/density diagnostic, is the topic 
of discussion for §5. The ratio of the two O VI features 
(A1038/A1032) is not provided in Table 6 but it is in all cases 
very close to 1:2. This is the value expected for an effectively- 
thin plasma where the flux ratio is proportional to the ratio of 
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the collision strengths for the two transitions, which for O VI is 
very close to 1:2 (Zhang, Sampson, & Fontes 1990). 

Plasma model fits to extensive X-ray observations of the 
three youngest stars in the sample (EK Dra, ir l UMa, and n 1 
Cet) indicate the presence of high-temperature coronal mate- 
rial (Giidel et al. 1997). As discussed by Feldman & Doschek 
(1991), Young et al. (2001) and Redfield et al. (2003), coro- 
nal emission lines of highly-ionized ions arising from forbid- 
den transitions are expected in the FUSE wavelength range. A 
coronal Fe XVIII forbidden line at 974.85 A was identified for 
the first time in FUSE spectra of Capella (Young et al. 2001). 
FUSE spectra of EK Dra, ir l UMa and k 1 Cet also clearly show 
this emission feature (see Figure 2 for an example). Very re- 
cently, Redfield et al. (2002) have reported the detection of a 
coronal Fe XIX Al 1 18 emission line in the FUSE spectra of sev- 
eral nearby cool stars. The FUSE spectra of EK Dra, tt 1 UMa 
and k 1 Cet show weak Fe XIX Al 1 18 features barely above the 
noise level. As pointed out by Redfield et al. (2003), this line 
can be corrupted by a nearby C I multiplet. Because of the poor 
S/N of the flux measurements for Fe XIX Al 1 18 and concerns 
with nearby line contamination have focused only on the clean 
Fe XVIII A975 line as a tracer of the hot coronal plasma. Inte- 
grated flux measurements have been carried out using the same 
method described above and are reported in Table 2. No trace of 
this feature is detected for the older stars in the sample down to 
the noise level of the spectra, and flux upper limits are therefore 
provided in Table 2. 

In addition to integrating the fluxes, we also obtained line 
centroid positions and FWHM measurements for the selected 
features in Table 2. The motivation for such measurements was 
the search for possible line shifts or broadenings that could arise 
from plasma outflows or extended material. A critical point 
when measuring line shifts is an accurate wavelength scale. 
Aperture centering problems caused by misalignments may re- 
sult in systematic velocity shifts among the different channels. 
Redfield et al. (2002) used ISM features to derive the abso- 
lute wavelength scale of their spectra. In our case ISM absorp- 
tion features are not measurable and we followed a different 
approach. Namely, we aligned the O VI AA1032,1038 doublet, 
which is available in all four channels, using the LiFl channel 
as the reference (this is the channel used for guiding). This pro- 
cedure should ensure relative velocities accurate to 5-10 km s -1 
in all channels. The velocity corrections were generally quite 
relevant with values up to 40 km s" 1 . With a corrected wave- 
length scale, we measured the line centroids and velocity shifts 
for several features. The results are given in Table 3. No signif- 
icant line shifts with respect to the photospheric velocity were 
observed down to the precision of the measurements. The ex- 
pected ^10 km s" 1 redshift of TR features (e.g., Redfield et al. 
2002) is only observed in the most active star of the sample, EK 
Dra. The rest of the stars actually exhibit slightly blueshifted 
TR features, albeit with very low statistical significance. 

FWHM measurements for all the emission lines in Table 
2 yielded relatively small widths of 0.20-0.25 A, which are 
equivalent to velocities of 60-70 km s" 1 . Only EK Dra, the 
most active of our targets, was found to have broader line pro- 
files. In this case, the measured FWHM was about 0.35 A 
or 110 km s _1 . (The actual rotational velocity of this star is 
v rot w 17 km s -1 .) A closer inspection of the C III A977 and O 
VI A 1032 line profiles revealed extended line wings in some of 
the targets. This is a well-known phenomenon (see, e.g., Linksy 
& Wood 1994) that is usually modeled by assuming double 
Gaussians with components of different widths. For the stars 



in our sample, only EK Dra shows unambiguous evidence of 
a broad component, whereas ir l UMa and k 1 Cet show some- 
what extended line wings but the two-component fits yielded 
inconclusive results. Two-Gaussian fits to the O VI A1032 line 
of EK Dra indicate a ratio between the broad component flux 
and the total flux of 0.38. Thus, extended line profiles ap- 
pear to correlate with stellar activity (and age), in agreement 
with the results of Redfield et al. (2002) and Wood, Linsky, & 
Ayres (1997). The presence of broad and narrow components in 
high-temperature lines is still not well understood. A plausible 
model (Wood et al. 1997) suggests that these may arise from 
high-velocity nonthermal motions during magnetic reconnec- 
tion events (microflares). 

3.1. Interstellar absorption corrections 

Before analyzing and modeling the line integrated fluxes, the 
possible effects of interstellar medium (ISM) absorption must 
be investigated. The C III A977 line and the H I Lyman se- 
ries lines are the most sensitive to attenuation effects because 
of possible superimposed ISM absorption. As shown in Table 
1, the stars in the sample are nearby and lie within 7-34 pc 
of the Sun and have negligibly small values of E(B-V) from 
ISM dust. However, even for the small ISM column densi- 
ties at these close distances (Nh ~ 10 18 cm" 2 ; Redfield & Lin- 
sky 2000), one expects some absorption in the lines mentioned 
above because they are among the strongest ISM features in the 
entire spectrum. Unfortunately, both the resolution and S/N of 
our FUSE spectra are not sufficiently high for reliable direct 
estimations of the (expected weak) ISM FUV absorptions. We 
therefore adopted an alternative procedure and carried out in- 
direct determinations of the flux corrections for the individual 
spectral features. 

To estimate the ISM column densities in the line of sight 
of our targets we made use of the local ISM observations and 
model by Redfield & Linsky (2002). These authors analyzed 
high-resolution observations of Fe, Mg and Ca features to de- 
termine column densities of these elements toward stars within 
100 pc of the Sun. One of our targets, k 1 Cet, is included in 
Redfield & Linsky's study and thus direct ISM measurements 
are available. For the other five targets we adopted the local 
ISM characteristics of neighboring stars with direct measure- 
ments: DK UMa for EK Dra and tt 1 UMa, HZ 43 for (3 Com, 
£ Dor for (3 Hyi, and 6 Cyg for 16 Cyg A. The C abundance 
was computed from the observed column densities of Mg II and 
Ca II and local ISM log abundances (relative to H) of -3.66, 
-5.58 and -7.64 for C, Mg II, and Ca II, respectively, from 
Wood et al. (2002b) and Redfield & Linsky (2000). The C III 
column density was subsequently obtained by adopting the lo- 
cal ISM ionization ratio of C Ill/C II«0.02 from Wood et al. 
(2002b). (Note that C II is the dominant ionization species in 
the ISM.) 

The ISM absorption equivalent widths were estimated from 
the C III column densities and the transition oscillator strengths 
through the curve of growth method (see Spitzer 1978). For 
the C III A977 transition, despite its strength, the ISM absorp- 
tion features for all our targets were found to be relatively weak 
(unsaturated) because of the small column densities. Then, us- 
ing the ISM component velocities and the observed C III line 
profile, we calculated the fraction of stellar flux absorbed by 
the ISM. Our results indicate that the flux corrections for the 
C III A977 line are very small (3-7%) for EK Dra, tt 1 UMa, k 1 
Cet, and (3 Com, and somewhat larger (1 1-13%) for (3 Hyi and 
16 Cyg A. 
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3.2. The H I Lyman series and total FUSE fluxes 

In addition to the features mentioned above, we attempted an 
estimation of the total integrated fluxes within the FUSE band- 
pass (920-1180 A). However, even when restricting to night- 
time spectra, contamination of the stellar H I Lyman features 
by geocoronal emission is very significant. This is because 
FUSE observations in this program were obtained through the 
large aperture to ensure, where possible, that the important C III 
977 A emission line was detected in the SiC channels. As a 
consequence, the Lyman geocoronal emission is much broader 
than for the medium aperture, reaching about ±100 km s" 1 near 
Ly/3. This emission overlies and in most cases completely dom- 
inates the stellar emission profile. The stellar signal is expected 
to be a self -reversed emission feature with a dark ISM absorp- 
tion core. 

To estimate the stellar fluxes for Ly/3 we took the difference 
between the night and day spectra, and assumed that the shape 
of geocoronal emission remains constant. A fraction of this 
geocoronal emission was then subtracted from the night spec- 
tra so that the residual flux at the location of the ISM feature is 
consistent with the flux predicted by a simulation of the stellar 
profile attenuated by the ISM. If emission peaks remained in 
the subtracted spectrum we assummed that the remaining sig- 
nal was stellar. Only k 1 Cet and n 1 UMa had any significant 
signal to satisfy these requirements. To estimate the shape of 
the underlying stellar emission profile, required to estimate the 
total stellar flux, we used the functional form of the "univer- 
sal" profile for strong partially coherent scattered lines derived 
by Gayley (2002). This allowed us to use the observed flux 
in the wings, which is relatively uncontaminated by geocoro- 
nal emission, to extrapolate towards line center and provide a 
total flux estimate. The model stellar profile was then attenu- 
ated by an ISM model (described above) and convolved with a 
PSF of R ~ 15000- 18000 to mimic the FUSE spectral reso- 
lution. From the fits we obtained total Ly/3 fluxes of 4.0T0" 13 
erg s" 1 cm" 2 and 2.4- 10~ 13 erg s" 1 cm" 2 for ir l UMa and k 1 Cet, 
respectively. Given the uncertainty of the spectra, the adopted 
ISM model and the assumed stellar profile shapes, we expect 
that the accuracy of estimated stellar flux is no better than a 
factor of 1 .5-2. Future observations of H Lya and MDRS ob- 
servations of Ly/3 are warranted and required to improve the 
quantification of the total stellar FUV emission. 

A correction to account for the rest of the H I Lyman series 
was derived by comparison with SOHO SUMER spectra for the 
Sun. We roughly estimate that the total flux contributed by the 
H I Lyman lines within the FUSE bandpass is about 1 .6 times 
the Ly f3 flux. Then, the total flux in the range 920-1180 A 
was calculated by adding the contributions from the individual 
features discussed above and the integrated H I Lyman series. 
We obtained values of 10.5-10" 13 erg s" 1 cm" 2 and 8.7-10" 13 
erg s" 1 cm" 2 for tt 1 UMa and k 1 Cet, respectively. Our inte- 
grations indicate that the H I Lyman features in this wavelength 
range (Ly/3, Ly7, etc) are important contributors (^40-60%) to 
the total flux of both ir 1 UMa and k 1 Cet. This is also presum- 
ably the case of the remaining solar-type stars in the sample but 
the weak stellar signal and strong geocoronal contamination of 
the FUSE observations prevented us from carrying out reliable 
measurements. 

4. EMISSION LINE FLUXES IN THE FUV 

One straightforward application of the measured line fluxes 
presented in Table 2 is the analysis of their evolution along the 



main sequence of a one solar mass star. To compare the abso- 
lute irradiances, we scaled the observed fluxes to surface fluxes 
at a radius of 1 R Q . To estimate the actual radius of each target 
we made use of the stellar luminosity (computed from the Hip- 
parcos distances and apparent magnitudes) and effective tem- 
perature (obtained from spectroscopic measurements in the lit- 
erature). The radii of the targets are included in Table 1 and 
using these we calculated the equivalent flux for a star of 1 R Q 
in radius. As can be seen in Table 1, the radii of all stars but 
(3 Hyi and 16 Cyg A are very close to 1 R Q (within 5-8%) and 
the correction was typically very small. The situation is some- 
what different for (3 Hyi because of its advanced evolutionary 
stage and its mass being about ~1.1 M . In this case the flux 
correction was quite significant. 

The adopted final stellar surface fluxes for the four emission 
features studied are provided in Table 4. Also in this table are 
the total integrated FUV fluxes in the FUSE wavelength range 
(920-1 180 A) for two of the targets (tt 1 UMa and k 1 Cet). All 
fluxes have been corrected for ISM absorption (see §3.1 and 
3.2). Before proceeding with a numerical analysis of these data, 
an illustrative visual impression of the surface flux changes for 
the targets stars can be seen in Figure 3. We have plotted de- 
tails of the spectra (flux corrected to the stellar surface) in a 
narrow 10- A wavelength interval that contains the O VI doublet 
at AA1032,1038. A rough comparison of the line variations as 
a function of stellar age does reveal a clear and definite trend or 
decreasing strength with increasing age. 

For a more quantitative analysis we used the integrated sur- 
face fluxes in Table 4, which are plotted in Figure 4 as a func- 
tion of the stellar rotation period. The features represented are 
C II A1037, C in A977, O vi A1032, and Fe xvin A975. The 
typical peak formation temperatures of the C II, C III, O VI, 
and Fe XVIII ions are approximately 20,000 K, 60,000 K, 
300,000 K, and 6 MK, respectively (Arnaud & Rothenflug 
1985). Thus, these lines probe a wide interval of plasma tem- 
peratures. As can be seen, the newly-detected coronal line of 
Fe XVIII is very important to our investigation because it ex- 
tends the coverage to the hot plasma component. Note the much 
smaller surface fluxes of (3 Hyi with respect to the younger so- 
lar analogs for all the features studied. Roughly, (3 Hyi's flux 
levels are comparable to, or slightly smaller than, those of to- 
day's Sun. In addition, Figure 4 also depicts the total integrated 
FUV fluxes in the FUSE bandpass for tt 1 UMa and k 1 Cet. 

The random measurement errors affecting the fluxes in Fig- 
ure 4 are discussed in §3. However, additional systematic un- 
certainties can be caused by the intrinsic stellar activity cycle. 
From an analogy with the Sun, we expect peak-to-peak cycle 
amplitudes in the FUV of about 30-50%. Since the target stars 
are observed at different stages of their activity cycles, these 
activity variations can be a source of scatter in the irradiance 
plot. Fortunately, this activity-related scatter is rendered negli- 
gible by the large relative differences between the fluxes of the 
targets (with factors of over 30 in flux). 

The fluxes of Figure 4 can be fit to good accuracy with power 
laws of different slopes. These are listed in Table 5, which con- 
tains the power-law slopes using both rotation period (P ro t) and 
age as independent variable. Note that only two FUSE mea- 
surements are available for the total flux within 920-1180 A. 
To constrain the fit we used the flux values for the quiet sun 
obtained through integration of the solar reference spectrum by 
Heroux & Hinteregger (1978) in the FUSE wavelength window. 
The resulting surface flux was 2.4- 10 4 erg s" 1 cm" 2 . 

Interestingly, the slopes not only change from line to line, but 
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also they show a clear trend. The power law becomes steeper 
as we move from cooler to hotter material. The most extreme 
trend is that of the flux for the coronal ion Fe XVIII. The re- 
sults indicate that the flux is reduced by a factor of over 1000 
when the rotation period increases ten-fold. To serve as a refer- 
ence, plasma of similar (or maybe slightly lower) temperature 
is probed when analyzing fluxes in the X-ray domain. Giidel et 
al. (1997) studied X-ray fluxes (0.1-2.4 keV) of about a dozen 
solar analogs covering also a wide period range. A power law 
fit to the data yielded a slope of -2.64 ±0.12, which is in very 
good agreement with our value for the decrease of the Fe XVIII 
flux. 

The results indicate that the emission from hotter plasma de- 
creases more rapidly than that from cooler plasma. Analogous 
behavior was observed by Giidel et al. (1997) when comparing 
the emission measures derived from observations in hard and 
soft X-rays. A possible scenario to explain the steeper flux de- 
crease for hotter plasma could be the weakening of the stellar 
dynamo as the stars spin down with age and the corresponding 
decline in the efficiency of heating mechanisms (e.g. flares) and 
in the strength of magnetic fields that confine the gas. 

It should be noted that our program stars, except for age and 
rotation period, have similar properties to the Sun (mass, radius, 
effective temperature, metal abundance). Importantly, these 
stars should all have similar convection zone depths, which to- 
gether with rotation, is an important parameter in most modern 
magnetic dynamo theories. Thus, the stars in our sample can 
serve as laboratories for testing the generation (and dissipation) 
of magnetic dynamo-related energy and heating in solar-type 
stars where rotation (angular velocity) is the only important 
variable. The observed power law dependencies on rotation of 
the various emission line fluxes should be important as inputs 
for testing and constraining stellar/solar dynamo theories. 

5. ELECTRON PRESSURES IN THE TRANSITION REGION 

The C III emission line at 977 A and the C III multiplet at 
1 176 A occur in the FUSE wavelength region. The ratio of the 
C III Al 176/A977 emission line fluxes has long been recognized 
to be a sensitive diagnostic for measuring the electron pressure 
(P e ) of hot (^60,000 K), optically-thin plasmas (e.g., Dupree, 
Foukal, & Jordan 1976; Keenan & Berrington 1985). FUSE of- 
fers an excellent opportunity to exploit the full potential of this 
important diagnostic (1Z = F^ m / F^n) because both C III fea- 
tures can be measured simultaneously. Furthermore, unlike the 
UV C III] A 1909 line, these FUV features are not contaminated 
by the stellar continuum or nearby line emissions. Details of the 
C III A977 and Al 176 emission features for two of the targets 
in the sample are shown in Figure 5. 

To compute 1Z ratios we adopted a simple C III model which 
includes the first 10 fine-structure levels. The energy levels 
are taken from Moore (1970), the highest level in the model 
being 2p 2 l S at 18219.88cm" 1 . The electric dipole oscilla- 
tor strengths are taken from the compilation of Allard et al. 
(1990), except for transitions calculated by Tachiev & Froese 
Fischer (1999; 2002, priv. comm.) and Fleming, Hibbert, & 
Stafford (1984). The electric quadrupole and magnetic dipole 
oscillator strengths are from Nussbaumer & Storey (1978) and 
Tachiev & Froese Fischer (1999) and the 2s2p 3 P| - 2p 2 1 Sq 
and 2s 2 l So~2s2p 3 P2 magnetic quadrupole strengths are from 
Shorer & Lin (1977) and Tachiev & Froese Fischer (1999), 
respectively. The electron collision strengths are taken from 
Berrington (1985) and Berrington et al. (1985). The proton 



collision rates for transitions within 2s2p 3 P° are from Ryans 
et al. (1998) and within 2p 23 P are from Doyle, Kingston, & 
Reid (1980). We solved the statistical equilibrium equations, 
under the assumption that the plasma is in steady state, with 
no photo-excitation from external sources. The ratios presented 
here are computed under the assumption that all the C III tran- 
sitions are optically thin. We have also computed the ratios 
under optically thick conditions using escape probabilities with 
hydrogen columns of 10 20 cm" 2 and assuming solar abundances 
(Grevesse & Sauval 1998). We find that 1Z is quite insensitive to 
optical depth effects under these conditions, in agreement with 
the findings of Bhatia & Kastner (1992). 

The A1176/A977 flux ratio of Be-like C III is well known 
to be sensitive to both electron density and temperature and is 
therefore sensitive to the shape of the emission measure dis- 
tribution between logr e = 4.5 and 4.9. The density sensitiv- 
ity arises from collisional thermalization of the lower levels of 
the 1 176 multiplet (2s2p 3 P°) which connect to the ground state 
(2s 2 1 S) by an intercombination electric dipole and a weak mag- 
netic quadrupole transition at 1908.7 A and 1906.7 A, respec- 
tively. The contribution functions for 1 176 A and 977 A contain 
the product of the collisional excitation which increases with T e , 
the ionization balance of C III which peaks near \ogT e = 4.85 
(Arnaud & Rothenflug 1985), and the shape of the emission 
measure distribution between log T e = 4.5 and 4.9 which is a 
strong declining function of T e , e.g., Jordan et al. (1987). The 
peak contribution is expected to lie below \ogT e = 4.85. For 
example for the Sun, Judge et al. (1995) find a peak as low 
as log T e = 4.5 for the 977 '-A transition, while MacPherson & 
Jordan (1999) find \ogT e = 4.7. For the solar analog transition 
regions, we consider that the C III emission lines are formed 
at a constant electron pressure, and we derive the ratio 1Z as 
a function of electron pressure rather than at a (less realistic) 
single temperature. 

The emission measure distribution at C III temperatures may 
change as a function of stellar magnetic activity. However, as 
pointed out by Jordan (2000), both the analyses of main se- 
quence stars by Jordan et al. (1987) and the similarity of the 
C IV/C II UV emission line ratios for a wide range of active 
coronal stars (Oranje 1986), indicate that the shape of the emis- 
sion measure distribution in the C III forming region approxi- 
mately follows the shape of the inverse of total radiative loss 
curve between logr c = 4.3 and 5.3. If the density sensitiv- 
ity of the ionization balance and abundance gradients can be 
neglected, then the relative electron pressures derived for our 
sample should be reliable. 

To quantify the effects of assuming differential emission 
measure distributions (DEM) with different slopes, we adopt 
a form 

DEM (T e ) = n e n H —^— cx T e ~ a (1) 
d\nT e 

where n e and n H are the electron and hydrogen densities, re- 
spectively, and r is the radial coordinate. Here we assume that 
the shape of DEM is similar to the emission measure over the 
C III 977 A formation region. Note that a ~ 2 is a reasonable 
approximation to the radiative power-loss function for this tem- 
perature range. For each pressure we integrate over temperature 
from log T e = 4.3 to 5.3, where the fractional abundance of C III 
exceeds 0.01, for each DEM to derive the total fluxes and de- 
rive 1Z. In Figure 6 we show the ratio 1Z for a=1.5, 2.0, and 2.5. 
This figure shows that the resulting ratio is not too sensitive to 
uncertainties in the shape of the emission measure distribution. 
Indeed, uncertainties in the collision strengths lead to greater 
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uncertainties in the electron pressure. 

The empirical C III ratios 1Z computed for the six stars in 
the sample are given in Table 6. Also provided are the TZ val- 
ues corrected for ISM absorption. From these one can estimate 
the electron pressures through the relationship discussed above. 
When doing so, we obtain P e values between 10 14 and 10 16 
cm -3 K, approximately. The actual electron pressures com- 
puted for each of the stars in the sample are given in Table 6. 
Both the TZ values and the derived plasma electron pressures 
are plotted in in Figure 7 as a function of the stellar rotation 
period. For comparison, the value for the quiet Sun was deter- 
mined to be 7^=0.29 (Dupree, Foukal, & Jordan 1976), and the 
active Sun (measured in solar active regions such as sunspots) 
has a C III ratio of 7Z=0.44 (Noyes et al. 1985; Doyle et al. 
1985). Note that the 1Z values determined for our stars corre- 
spond to the integrated stellar disk average, thus including both 
active and inactive regions, weighted towards region of large 
areal filling factor and high n e . Interestingly, the observed data 
for the target stars suggest a strong correlation between the ac- 
tivity level of the star (= age or P rot ) and 1Z or P e . 

The elusive density of the TR has been probed for the first 
time along the evolutionary path of a solar-type star, from 
ZAMS to TAMS. Our results show that the electron pressure 
of the TR decreases monotonically as the star spins down (and 
magnetic dynamo activity decreases). From the fit in Figure 7 
we infer a power-law relationship with a slope of ~ -1.7. This 
suggests that the electron pressure of the solar TR may have 
decreased by about a factor of ^40 since the beginning of its 
main sequence evolution. A likely explanation for the greater 
plasma electron pressure for the young and active (rapidly ro- 
tating) stars is a stronger magnetic confinement of the emitting 
material. 

Most notably, the power-law slope for the electron pressure 
and the C III 977 A flux are very similar, namely ~ —1 .65. Note, 
however, that the flux power-law relationship has been com- 
puted from the integrated flux and it thus represents an average 
over the stellar surface. If, however, the emitting TR plasma 
only occupies a fraction (A; A < 1) of the surface, known as 
filling factor, then the specific flux emitted by the active regions 
will be correspondingly larger. The pressure power-law does 
not depend on the area factor if only one atmospheric compo- 
nent dominates. Thus, assuming that the same fraction of C III 
977 A is emitted over the same d\\\T e interval for our sample 
then we can adopt 



F A=i 2 dr 

oc / 

A e d\\\T e 



(2) 



where F A= \ represents the surface integrated flux. Using the 
power-law relationships derived above, one can write 



d\nT e 
dr 



cxAf£ t 77 



(3) 



for our solar proxies 

log A ~ -0.86 R ^ -0.068 P mt 
and thus, the temperature gradient can be expressed as 



d\nT e 
dr 



Kj p-1.77 10 -0.068P ro 



(4) 



(5) 



as a function of the rotation period. Thus, with empirically 
measured filling factors for the solar proxies or with a reli- 
able relationship between A and P rot (or Rossby number 5 ) the 
relationship above can yield an estimate of the dependence of 
the electron temperature gradient with age. For example, using 
the empirical fit of Montesinos & Jordan (1993) that relates the 
magnetic filling factor with the Rossby number, we can write 

5 The Rossby number (i? ), which is proportional to the rotation period of the 
zone, is commonly used as a measure of the activity level of a star. For our targets, 
Rossby number is simply proportional to P rot . 



A representation of the relative variations of the temperature 
gradient, the electron pressure and the specific C III A977 flux 
(i.e., flux divided by filling factor) are provided in Figure 8. 
Note that these relationships are only represented for rotation 
periods shorter than that of the Sun, where the filling factor is 
better defined. Strikingly, the specific C III A977 flux remains 
relatively constant with rotation period: as the filling factor in- 
creases, the total stellar flux increases correspondingly. How- 
ever, since the electron density increases with decreasing ro- 
tation period, the temperature gradient must also increase to 
keep the specific flux relatively constant. It is worth noting 
that, given the large change in the implied temperature gradi- 
ent across the sample, one or more of the several assumptions 
on which Eq. 3 is based may no longer hold. 

6. FUV IRRADIANCES AND EFFECTS ON PALEO-PLANETARY 
ATMOSPHERES 

The relationships obtained from the solar proxies indicate 
that the total solar FUSE/FUV flux was about twice the present 
value 2.5 Gyr ago and about 4 times the present value about 3.5 
Gyr ago. Note also that the FUSE/FUV flux of the Zero- Age 
Main Sequence Sun could have been stronger than today by as 
much as 50 times. 

As discussed by Canuto et al. (1982, 1983), Luhmann & 
Bauer (1992), Ayres (1997), Guinan, Ribas & Harper (2002) 
and others, the strong FUV and UV emissions of the young, 
more active Sun could have played a major role in the early 
development and evolution of planetary atmospheres - espe- 
cially those of the terrestrial planets. The expected strong X- 
ray-UV irradiance of the young Sun can strongly influence the 
photochemistry and photoionization (and possible erosion) of 
the early planetary atmospheres and even surfaces (in the case 
of Mercury, Moon and Mars) and also may play a role in the 
origin and development of life on Earth as well as possibly on 
Mars. For example, Canuto et al. (1982, 1983) discuss the 
photochemistry of O2, O3, CO2, H2O, etc, in the presumed out- 
gassed C02-rich early atmosphere of the Earth. Ayres (1997) 
discusses the effect of the young Sun's increased ionizing X- 
ray-UV flux, and possible accompanying enhanced solar wind, 
on the erosion of the early atmosphere of Mars about 3-4 Gyr 
ago. Also, Lammer et al. (2003) utilize irradiance data from the 
"Sun in Time" project and solar-wind estimates by Wood et al. 
(2002a) to evaluate the mechanisms for loss of water from Mars 
and study the implications for the oxidation of the Martian soil. 

Similarly, our data can also provide insights into the so-called 
Faint Sun Paradox. The paradox arises from the fact that stan- 
dard stellar evolutionary models show that the Zero- Age Main 
Sequence Sun had a luminosity of ^70% of the present Sun. 
This should have led to a much cooler Earth in the past while 
geological and fossil evidence indicate otherwise. A solution 
to the Faint Sun Paradox proposed by Sagan & Mullen (1972) 

star and inversely proportional to the turnover time at the base of the convective 
which have very similar masses and spectral types, one can safely assume that the 
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was an increase of the greenhouse effect for the early Earth. 
The gases that have been suggested to account for this enhanced 
greenhouse effect are carbon dioxide (CO2), ammonia (NH3) or 
methane (CH 4 ). However, recent results for atmospheric com- 
position of the early Earth (Rye, Kuo, & Holland 1995) are in 
conflict with the high levels of CO2 and H2O needed to ex- 
plain the stronger greenhouse effect. Moreover, ammonia is a 
likely candidate except that it is quickly and irreversibly pho- 
todissociated. However, Sagan & Chyba (1997) proposed that 
hydrocarbon-based aerosols could shield the ammonia from 
damaging incoming solar radiation. Another alternative expla- 
nation for the enhanced greenhouse effect is from atmospheric 
methane as discussed by Pavlov et al. (2000). Methane is a 
very strong greenhouse gas (10 3 — 10 4 times stronger than CO2) 
and the necessary amount of CH4 could have been provided by 
methanogenic bacteria in the proposed early Earth's reducing 
atmosphere. 

In addition to the stronger greenhouse effect, the expected 
higher X-ray-UV irradiance of the young Sun should heat up 
the upper atmosphere of the early Earth. Although the stronger 
high-energy solar radiation cannot by itself explain the Faint 
Sun Paradox, the photoionization and photodissociation reac- 
tions triggered could play a major role in what greenhouse gases 
are available. For example, the high levels of FUV-UV radia- 
tion of the young Sun could strongly influence the abundances 
of ammonia and methane in the pre-biotic and Archean plane- 
tary atmosphere some 2-4 Gyr ago. Similarly, the photochem- 
istry and abundance of ozone (O3) is of great importance to 
study life genesis in the Earth. Ozone is an efficient screening 
mechanism for the enhanced UV radiation of the young Sun, 
thus protecting the emerging life on the Earth's surface. 

To fully evaluate the influence of FUV radiation in the de- 
veloping planetary system one must first estimate the flux con- 
tribution of the strong H I Lya feature, which is not included 
in the FUSE wavelength range. Preliminary estimates using 
HST/STIS spectra indicate that Lya A1216 contributes a sig- 
nificant fraction (up to 90%) of the total FUV flux. Further 
observations and measurements in the near future will allow us 
to obtain this missing piece of information and complete the 
FUV irradiance study. 

7. CONCLUSIONS 

In this study we have utilized FUV spectra acquired with the 
FUSE satellite to investigate the FUV emission characteristics 
of solar analogs. Our six targets are well-known G0-5 solar- 
type stars especially selected to serve as proxies for the Sun at 
different ages, nearly covering the entire solar main sequence 
lifetime from 130 Myr to ~9 Gyr. Here we have focused on 
two aspects of the "Sun in Time" program which are the study 
of TR plasma electron pressures (using C III Al 176/A977 line 
ratio diagnostics) for stars that differ only in rotation period 
(and age) and the evolution of irradiances for specific features 
covering emitting plasma with temperatures from ~ 10 4 K to 
- 10 7 K. 

To analyze the plasma density of the TR we have used a the- 
oretical relationship between the ratio of C III Al 176 and A977 
emission line fluxes and the electron pressure (P e ) of the mate- 
rial responsible for the emission. Fortunately, both C III tran- 
sitions are within the FUSE wavelength range so that they can 
be measured simultaneously. Our results indicate a power-law 
relationship between the electron pressure and the stellar rota- 
tion period and overall magnetic -related activity (both related 
to age). The slope of this relationship has been found to be 



~— 1.7, which suggests that the electron pressure of -10 5 -K 
material in the Sun has decreased significantly since the begin- 
ning of its main sequence evolution. The higher values of elec- 
tron pressure found for the more active stars are best explained 
by stronger magnetic confinment of the TR plasma. 

The measured fluxes for four emission features - C II A 1037, 
C III A977, O VI A 1032, and Fe xvm A975 - were referred to a 
radius of 1 Rq and also corrected for ISM absorption whenever 
necessary. The typical formation temperatures of the studied 
features are 20,000 K, 60,000 K, 300,000 K, and 6 MK, re- 
spectively. Our analysis indicates that the evolution of these 
fluxes with stellar rotation period or age can be accurately fit 
with power-law relationships of different slopes. Interestingly, 
the slopes not only change from line to line, but also they show 
a clear trend: The power law becomes steeper as we move from 
cooler to hotter plasmas, with the most extreme trend being that 
of the coronal ion Fe XVIII. Also, the evolution total integrated 
flux in the FUSE wavelength range (920-1180 A) can be fit 
with a power-law relationship with a slope of ~— 1.7, which 
indicates a factor of ~50 flux decrease along the solar main 
sequence evolution as the magnetic dynamo activity decreases. 

The high levels of FUV line emission fluxes (and related 
high-energy emission) of the early Sun could have played a 
crucial role in the photochemistry and photoionization of ter- 
restrial and planetary atmospheres. To address this point we are 
completing spectral irradiance tables covering 1 A to 3200 A 
for our program stars that can serve as input data for evolution 
and structure models of the paleo-atmospheres of the Solar Sys- 
tem planets. The FUSE observations fill a critical wavelength 
and energy gap in the "Sun in Time" program and complement 
observations of the same stars in the X-ray and EUV regions 
(corona) made with ROS AT, SAX, ASCA, XMM, Chandra, and 
EUVE, and in the UV (TR and chromosphere) made using IUE 
and HST The only missing piece of information is the evo- 
lution of the irradiance of the strong chromospheric H I Lya 
A1216 FUV feature. We are currently carrying out this part of 
the study and we expect to complete it shortly. Full spectral 
irradiance tables for five solar proxies are thus nearing comple- 
tion and will be made available in a forthcoming publication 
(Ribas & Guinan 2003, in preparation). 

We anticipate that these irradiance results will be important 
for the study of paleo-atmospheres of the Solar System plan- 
ets. In particular, preliminary analyses indicate that the high 
X-ray and EUV emission fluxes of the early Sun could have 
produced significant heating of the planetary exospheres and 
upper-atmospheres thus enhancing processes such as thermal 
escape. The early Sun's strong FUV and UV fluxes penetrate 
further into the atmosphere and probably influenced the photo- 
chemistry of, e.g., methane and ammonia, which are important 
greenhouse gases. Thus, this study has strong implications for 
the evolution of the pre-biotic and Archean atmosphere of the 
Earth as well as for the early development of life on Earth and 
possibly on Mars. 

We thank Seth Redfield for making available to us the data 
from the Colorado LISM model, which we have used to cor- 
rect the measured fluxes for ISM absorption. The referee, Tom 
Ake (Johns Hopkins Univ.), is thanked for helpful comments 
and suggestions that led to significant improvements. We ac- 
knowledge with gratitude the support for the "Sun in Time" pro- 
gram from NASA-FUSE grants NAG 5-08985, NAG 5-10387, 
NAG 5-12125 and also from NSF-RUI grant AST-00-71260. 
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Table 1 

Relevant data for the "Sun in Time" FUSE targets and the Sun. 



Name 


HD 


Sp. 


d 


?eff 


Mass 


Radius 


Prot 


Age 


Age 






Typ. 


(pc) 


(K) 


(Mq) 


(R©) 


(d) 


(Gyr) 


indicator 


EK Dra 


129333 


GOV 


33.9 


5818 


1.07 


0.95 


2.75 


0.13 


Pleiades str. 


tt 1 UMa 


72905 


G1.5 V 


14.3 


5840 


0.98 


0.97 


4.68 


0.3 


UMa str. 


k' Cet 


20630 


G5 V 


9.2 


5700 


1.01 


0.94 


9.2 


0.65 


Prot-Age rel. 


/3Com 


114710 


GOV 


9.2 


5950 


1.11 


1.08 


12.4 


1.6 


P ro i-Age rel. 


Sun 




G2 V 


1 AU 


5777 


1.00 


1.00 


25.4 


4.6 


Isotopic dating 


/3Hyi 


2151 


G2IV 


7.5 


5800 


1.09 


1.88 


~28 


6.7 


Isochrones 


16 Cyg A 


186408 


G2 V 


21.6 


5740 


0.99 


1.27 


~35 


^9 


Isochrones 
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Table 2 

Measured Integrated Fluxes for Selected Emission Features 



Name 


ObsID 


Date Obs. 


tcxp (ks) 


-ttCiii 
" r A977 


ttCiii 
^ A1176 


-'A1032 


-r-Ovi 
-^A1038 


A1037 


•r-Ftjxviii 
- r A975 


EKDra 


C1020501 


2002 May 14 


24.2 


11.1 ±0.7 


7.6±0.4 


7.3±0.3 


3.4±0.2 


0.47±0.07 


0.81±0.19 


tt 1 UMa 


B0780101 


2001 Dec 5 


16.5 


14.7±0.7 


9.3±0.5 


9.6±0.3 


4.8±0.2 


0.74±0.09 


0.65±0.26 


k 1 Cet 


A0830301 


2000 Sep 10 


13.0 


16.4±0.7 


10.7±0.4 


12.4±0.4 


6.2±0.2 


1.07±0.12 


0.44±0.14 


Com 


A0830401 


2001 Jan 26 


15.0 


11.0±0.5 


5.5±0.2 


5.9±0.2 


2.8±0.1 


0.63±0.08 


<0.06 b 


/3Hyi 


A0830101 


2000 Jul 1 


18.0 


15.1 ±0.5 


5.8±0.2 


6.1 ±0.2 


2.8±0.1 


1.02±0.09 


<0.04 b 


16 Cyg A 


C1020101 


2002 Jul 30 


35.0 


0.66±0.14 


0.20±0.09 


0.43±0.05 


0.20±0.03 


<0.05 b 


<0.04 b 



a Fluxes in units of 10 14 erg s 1 cm" 
b l-(T upper limits 



Table 3 

Radial Velocity Measurements for Selected Features in the FUSE Spectra 3 



Name 


Phot 


CO 


A 1037 


Cm 


A977 


OVI 


A 1032 


FeXVIIlATO 




Hel 


Hel 


Shift 


Hel 


Shift 


Hel 


Shift 


Hel 


Shift 


EK Dra 


-30.5 


-23 


8 


-15 


16 


-25 


6 


-26 


5 


tt 1 UMa 


-12.0 


-23 


-11 


-24 


-12 


-19 


-7 


-51: 


-39: 


k 1 Cet 


19.9 


12 


-8 


16 


-4 


16 


-4 


1 


-19 


/3 Com 


6.1 


3 


-3 





-6 


4 


-2 






/3Hyi 


22.7 


6 


-17 


20 


-3 


16 


-7 






16 Cyg A 


-25.6 


-40: 


-14: 


-29 


-3 


-22 


4 







a All velocities in units of km s" 1 . Phot=stellar photosphere; 
Hel=heliocentric velocity; Shift=Doppler shift with respect to stellar pho- 
tosphere. 



Table 4 

ISM- and Radius-corrected Surface Fluxes for the Five of the Targets with High S/N Data' 



Name 


/■Cn 
J A1037 


/•Cm 
J A977 


J A 1032 


/■Fexvin 
J A975 


/■FUV 
J 920-1 180 


EK Dra 


1.17±0.18 


28.4±1.7 


18.2±0.8 


2.02±0.51 




tt 1 UMa 


0.32±0.04 


6.52±0.33 


4.10±0.12 


0.28±0.11 


45 ±20 


k 1 Cet 


0.20±0.02 


3.30±0.15 


2.32±0.07 


0.082±0.025 


16.4±3.8 


Com 


0.090±0.011 


1.64±0.07 


0.842±0.030 


<0.0085 




/3Hyi 


0.032±0.003 


0.543±0.019 


0.190±0.007 


<0.0()12 




16 Cyg A 


<0.03 


0.423±0.090 


0.246±0.028 


<0.02 





a Surface fluxes in units of 10 4 erg s 1 cm 2 
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Table 5 

Flux decrease power law slopes for different wavelength regions and ions. 



Wav. region/ion Plasma T (K) Slope (P ro ,) a Slope (age) b 



FUV (920-1180 A) 


~ 10 4 -10 5 


-1.74 


-1.03 


Co 


2- 10 4 


-1.49 


-0.88 


C III 


6- 10 4 


-1.61 


-0.95 


VI 


3 ■ 10 5 


-1.82 


-1.07 


X-rays (0.1-2.4 keV) c 


2- 10 6 


-2.64 




Fe XVIII 


6-10 6 


-3.2 


-1.9 



a Power law slope corresponding to the following func- 
tional form: flux = flux (P TOt ) a , where a is the exponent 
listed in the column 

b Power law slope corresponding to the following func- 
tional form: flux = flux (age)* 3 , where (i is the exponent 
listed in the column 

c From Gudel et al. 1997 



Table 6 

C in Al 176/A977 flux ratio and electron pressure of the transition region 



Name 






EKDra 


0.68±0.05 


0.67±0.05 


7T 1 UMa 


0.63±0.05 


0.61±0.04 


k 1 Cet 


0.65±0.04 


0.60±0.04 


(3Com 


0.50±0.03 


0.48±0.03 


/3Hyi 


().38±0.02 


0.33±0.02 


16 Cyg A 


0.30±0.15 


0.27±0.13 



logP„(cm- 3 K) 



>16.0 



15 4+ 06 
15 3+ 05 

;8:f 



14 

14 1 +(M 



a-T-> — ■77cm / -77cm 

' V — A 1 176/ A977 

b C III A977 flux corrected for ISM absorption 
(see text) 
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FIG. 1 . — SiC2A (top), LiFIA (middle) and LiF2A (bottom) night-time spectra of the target star k 1 Cet (see Table 1 for information on the stellar properties). The 
most prominent emission features are labeled with the name of the corresponding ion. 
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FIG. 2. — Detail of the FUSE spectrum of re 1 Cet in the wavelength region around the coronal Fe XVIII emission feature. Also present are the stellar C III emission 
line and some "airglow" features. 
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FIG. 3. — Comparison of the surface fluxes for all targets in our FUSE programs in the region around the O VI AA1032,1()38 doublet. The spectra have been 
zeropoint-shifted using integer multiples of 10 5 erg s -1 cnT 2 A -1 to avoid confusion. Note the obvious trend of decreasing flux with increasing stellar age. The 
feature located near A = 1039 A is of geocoronal origin. 
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FIG. 4. — ISM- and radius-corrected surface fluxes vs. rotation period for our six FUSE targets. Plotted here are flux measurements for four features corresponding 
to transitions of C II, C III, O VI and Fe XVIII, which trace plasmas with temperatures of 20. 000 K, 60, 000 K, 300. 000 K and 6 MK, respectively. Also included 
are measurements of the integrated FUV (920-1 180 A) flux in the FUSE wavelength range for two of the targets. In this case, the integrated solar flux in the same 
wavelength interval was used to constrain the fit (see text). The error bars of the measurements are plotted but smaller than the size of the symbols in most cases. 
The straight lines are power law fits to the measured data and the corresponding slope coefficients are printed at the right end. 
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FIG. 5. — Details of the C III A977 and Al 176 emission features in the k 1 Cet and j3 Hyi FUSE spectra. The C III Al 176 feature is a multiplet that comprises six 
transitions. The flux scale in the C III Al 176 panels has been enhanced by a factor of 4 for illustration purposes. 
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FIG. 6. — The flux ratio 1Z for constant pressure models with three differential emission measure distributions (DEM). In the C III formation region the DEM 
decreases with increasing T e . The gradient of the DEM increases from the top curve down; top (dotted line) a = 1.5, middle (solid line) a = 2.0 and bottom (dashed 
line) a = 2.5 - where the DEM is given by Eq. 1. The characteristic formation temperature decreases with increasing a and gives rise to the to a lower ratio at a 
given pressure. The difference in the curves is small compared to the uncertainties in the observed ratio. 
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FIG. 7. — Top: Observed values of the flux ratio 72. = C III Al 176/A977 as a function of the stellar rotation period (which is an age indicator) for the stars in the 
sample. The typical TZ values for the Sun (corresponding to quiet and active regions) are also plotted. Bottom: Electron pressures derived from the C III Al 176/A977 
diagnostics using the theoretical calculations in Figure 6. The straight line is a power law fit (all data excluding the Sun) and the resulting expression is explicitly 
given in the panel. 
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FIG. 8. — Variations in the temperature gradient, electron pressure and specific C III A977 flux scaled to the solar values. The derivation of the relationships plotted 
is discussed in §5. Note the steep decrease in the temperature gradient with increasing rotation period (or age) and the slowly varying specific C III A977 flux (i.e. 
flux in active regions computed through an empirical relationship between the filling factor and P ro t). The plots stop at the rotation period of the Sun because of 
increasing uncertainties resulting from the filling factor-P ro t relationship. 



